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INTRODUCTION 

In  low  Earth  orbit,  energetic  collisions  between  molecules  from  spacecraft  emissions  and  atmosphere  are  observed 
to  produce  a  luminescent  interaction  layer  [1,  2].  These  collisions  deposit  a  substantial  amount  of  energy  into  the 
internal  modes  of  the  molecules.  The  rate  at  which  this  energy  is  deposited  and  the  distribution  of  that  energy  among 
translational,  rotational,  and  vibrational  degrees  of  freedom  govern  many  properties  of  the  interaction  layer,  including 
its  temperature  distribution,  shape,  and  spectral  radiance. 

Direct  simulation  Monte  Carlo  (DSMC)  models  [3]  are  among  the  important  tools  for  predicting  the  interactions 
between  the  atmosphere  and  spacecraft  surfaces  or  gases  emitted  by  spacecraft  in  low  Earth  orbit  [4-6].  Recent 
experimentally-determined  vibrational  transition  chemistry  (such  as  for  O  +  CO  [7])  now  enables  use  of  an  internal 
energy  model  with  explicit  transition  cross-sections  between  ground  and  explicitly-enumerated  quantized  excited 
vibrational  states  which  we  will  explore  in  this  work.  The  SOCRATES1  and  SOCRATES-P2  codes  [8-10]  are  DSMC 
codes  that  are  specialized  for  computing  spacecraft  contamination  flow  fields  [11],  radiative  byproducts  of  space 
vehicle  effluents  and  exhausts  [12-14]  and  their  spectra,  and  for  use  in  analysis  of  molecular  beams  [15,  16].  The 
SOCRATES  code  was  originally  developed  in  the  1980s  and  uses  computationally  economic  approaches  such  as  a 
variable  hard-spheres  (VHS)  interaction  model  [17]  with  a  single  global  ft)  parameter,  a  global  constant  probability 
for  inelastic  collisions,  and  a  single  per-species  constant  for  the  number  of  internal  degrees  of  freedom.  The  parallel 
version  of  the  code  has  several  additional  models  for  total  cross  sections  and  energy  redistribution.  Despite  the  apparent 
simplicity  of  the  original  models  in  SOCRATES,  they  are  highly  successful  in  their  primary  applications.  In  this  work 
we  use  DSMC  calculations  in  a  simplified  geometry  to  explore  the  sensitivity  of  flow  results  to  the  internal  energy 
models  by  comparing  flows  computed  using  several  internal  energy  and  energy  transfer  schema,  including  a  model 
with  detailed  quantized  vibrational  transitions. 


MODELS 

Several  models  are  pertinent  to  the  transfer  of  internal  energy,  including  the  cross-section  model,  the  internal  energy 
models  in  each  species,  and  energy  redistribution  algorithms.  We  begin  with  the  different  internal  degrees  of  freedom 
among  which  energy  is  redistributed  during  a  collision. 

Translational  degrees  of  freedom  are  governed  by  the  cross-section  model.  For  this  work,  the  sole  potential  model  is 
the  variable  hard  spheres  (VHS)  model  [17],  that  allows  the  effective  cross-section  of  a  molecule  to  vary  with  collision 
velocity  through  a  power-law  relation. 

d  =  (/ref  (c,- ref / Cr) (<B_ * /2)  ( 1 ) 


1  Spacecraft/Orbiter  Contamination  Representation  Accounting  for  Transiently  Emitted  Species 

2  SOCRATES,  parallel  version 
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The  ft)  exponent  is  significant  to  internal  energy  redistribution  because  there  are  (5  —  2co)  translational  degrees  of 
freedom  in  a  collision  involving  two  VHS  molecules.  The  VHS  parameters  are  set  using  viscosity  data  [18-20] 
according  to  the  methods  of  Bird  [3],  sec.  3.5  and  4.3. 

The  next  category  of  degrees  of  freedom  is  for  molecular  internal  energy.  The  simplest  internal  energy  model  is  a 
lumped  constant  model.  In  this  model,  all  non-translational  internal  energy  (rotational,  vibrational,  and  electronic)  is 
lumped  into  C,  internal  degrees  of  freedom,  such  that  on  average,  £)  =  (£/2)k7]  and  'C  is  constant.  C  relates  directly  to 
heat  capacity  Cp. 

C=^-5  (2) 

This  is  not  to  suggest  that  Cp  data  could  be  used  outright  to  assign  C,  but  rather  that  Cp  gives  (Q  where  the  ensemble 
is  considered  is  one  in  equilibrium  with  the  temperature  of  the  tabulated  Cp.  Consider  the  CO  molecule  as  described 
in  the  high-temperature  thermodynamic  database  of  McBride  and  Gordon  [18].  At  low  temperature  CO  has  a  C,  of 
2.0,  but  approaching  ©v,  C  increases  to  4.0.  As  the  energy  regime  of  our  calculation  is  commensurately  high,  we  will 
examine  each  of  these  values. 

In  a  collision,  the  Larsen-Borgnakke  [21]  method  is  used  to  assign  energy  to  every  degree  of  freedom  in  the  collision 
complex.  The  method,  however,  says  nothing  about  the  mechanism  for  or  the  rate  at  which  these  processes  occur — 
only  the  outcome.  The  rate  is  governed  by  an  externally  assigned  conditional  probability  that,  given  a  collision  occurs, 
it  is  inelastic. 

This  conditional  probability,  or  fraction  of  collisions  that  are  inelastic  (ft)  has  in  the  past  been  set  arbitrarily  to 
50%  for  lack  of  better  information.  Higher-fidelity  treatments  for  the  vibrational  and  rotational  collision  number  [22- 
24]  may  be  beneficial  because  they  add  a  local  collisional  energy  dependence  to  the  collision  number.  The  focus  of 
this  study,  in  the  area  of  energy  transfer  rate,  is  limited  however  to  comparing  a  global  constant  fraction  of  inelastic 
collisions  with  the  chemical-reaction  based  vibrational  transfer  kinetics. 


Explicit  Models  for  Vibrational  Exchange  Kinetics 

When  explicit  information  regarding  a  certain  type  of  energy  transfer  is  available,  it  is  preferable  to  use  it  instead 
of  the  Larsen-Borgnakke  and  collision-number  methods  for  that  particular  degree  of  freedom.  One  way  to  model 
vibrational  energy  transfer  is  to  treat  vibrational  transitions  as  explicit  chemical  reactions.  Under  this  methodology, 
vibrational  degrees  of  freedom  are  treated  explicitly  and  separately.  For  example,  in  a  lumped  internal-energy  model 
of  CO  that  usually  has  C,  =  4.0,  it  would  instead  remain  with  C  =  2.0;  explicitly  excited  species  [CO(v)}  such  that 
v'  >  0  are  added  explicitly  through  a  chemical  reaction  model.  Now,  at  the  cost  of  complexity,  we  have  control  of  the 
energy  dependent  rates  of  the  vibrational  transitions. 

Chemical  reaction  rate  constants  are  generally  parameterized  in  the  Arrhenius  form. 

k  =  AT,,exp(Ea/kT)  (3) 

These  rate  constants  are  measured  in  equilibrium  gas,  with  a  distribution  of  collision  velocities.  For  use  in  DSMC,  the 
reaction  is  Laplace-inverse  transformed  [17]  to  yield  the  following  cross-section: 

VrOr  =  VS  (1  +  *ij)Ay/\-Ea/Ec(Ec  -  Ea)n  /(2k”  r(n  +  1 .5))  (4) 

where  vr  is  the  relative  velocity  of  the  collision  pair,  ar  is  the  reactive  cross-section,  8(  /  is  one  if  the  colliding  molecules 
are  the  same  species  and  zero  otherwise,  and  Ec  is  the  translational  collision  energy.  Given  that  a  pair  of  species  collide, 
the  probability  of  reaction  is  a,  /Ovhs- 

Two  mechanisms  are  considered:  a  two-state  mechanism  and  a  multi-state  mechanism.  Although  we  include  the 
radiative  chemistry  here,  we  will  turn  it  off  for  some  of  the  calculations  to  separate  out  its  influence.  The  two-state 
mechanism  involves  ground  and  excited  CO: 

CO  +  M  CO(*)  +M 
CO(*)  — >  CO  T  photon 

The  2-state  model  uses  an  excitation  cross-section  fitted  to  the  data  of  Braunstein  and  Duff  [25],  wherein  A  = 
2.81  x  1(T25  m3s_1K_",  n  =  2.2,  and  Ea  =  4.25  x  10  20  J.  For  the  reverse  reaction,  the  same  A  and  n  is  used,  but 
the  activation  energy  (Ea)  is  assigned  to  zero. 


TABLE  1.  Matrix  of  parameters  surveyed  for  the  system 
O+CO.  Off-nominal  parameters  are  bold-typeset  for  emphasis. 


Label 

CO 

C 

fi 

Explicit  Vibrations 

Radiative 

N 

0.75 

2.0 

0.5 

none 

no 

A 

0.75 

2.0 

0.0 

none 

no 

B 

0.75 

2.0 

0.1 

none 

no 

C 

0.75 

2.0 

1.0 

none 

no 

D 

0.75 

4.0 

0.5 

none 

no 

E 

0.90 

2.0 

0.5 

none 

no 

F 

0.75 

2.0 

0.5 

2-state 

no 

G 

0.75 

2.0 

0.5 

multistate 

no 

H 

0.75 

2.0 

0.5 

2-state 

yes 

J 

0.75 

2.0 

0.5 

multistate 

yes 

The  multistate  mechanism  involves  several  reactions: 

CO(v)+M  ±5  CO(v')+M 
CO(v'  =  v+l)  — >  CO  (v)+  photon 

where  in  the  excitation  reaction,  v'  may  be  much  higher  than  the  original  v  state.  The  preceding  mechanisms  are  treated 
as  independent  of  the  J  state.  For  the  multistate  model,  only  the  reactive  cross  sections  of  Figure  6a  in  Brunsvold  et  al. 
[7]  were  used.  Brunsvold  et  al.  separately  determined  the  cross  sections  for  vibrational  excitation  due  to  reaction  (O 
substitution)  and  to  inelastic  collisions.  Using  the  cross-sections  for  the  reactive  channels  only  was  a  simplification  for 
convenience,  based  on  the  finding  of  Brunsvold  et  al.  that  the  reactive  and  inelastic  vibrational  excitation  cross-sections 
were  comparable.  The  simplification  will  lead  to  underestimating  the  incidence  of  vibrational  excitations  perhaps  by  a 
factor  of  two. 

In  this  context,  the  relevant  difference  between  the  two-state  and  the  multi-state  mechanisms  is  that  in  the  two-state 
mechanism  limits  the  amount  of  energy  deposited  into  the  vibrational  mode  to  one  quantum,  whereas  the  multistate 
model  allows  multiple  quanta. 


CALCULATIONS 

Our  calculations  involve  a  highly-simplified  system  in  which  there  are  two  translationally-cold  flows  impinging  at  right 
angles.  This  is  analogous  to  a  crossed  molecular  beam  experiment  or  a  thruster  firing  at  90°  to  the  travel  direction. 
Figure  1  shows  the  flow  schematically.  All  four  boundary  conditions  are  vacuum  boundaries — molecules  that  cross  the 
boundary  are  removed  from  the  simulation.  Flow  is  introduced  from  the  left  and  bottom  boundaries.  Though  decidedly 
imperfect,  this  geometry  is  a  convenient  tool  for  qualitatively  surveying  the  sensitivity  of  the  scattering  behavior  to  the 
different  internal  energy  transfer  and  internal  energy  models.  In  the  preceding  section,  we  identified  several  parameters 
that  might  affect  the  shape  of  a  flow  feature:  CO,  C  fi  (the  fraction  of  collisions  that  are  inelastic),  and  explicit  vibrational 
kinetics.  We  survey  these  conditions  using  the  system  of  atomic  oxygen  (O)  and  carbon  monoxide  (CO)  at  a  center-of- 
mass  relative  velocity  of  8000  m/s.  Table  1  shows  the  conditions  explored  in  this  survey.  The  nominal  conditions,  with 
label  “N,”  are  for  ©  =  0.75,  C  =  2.0,  f,  =  0.5,  and  no  radiation  or  explicit  vibrational  energy  transfer.  This  condition 
best  corresponds  to  a  low-energy  CO  collision. 

All  calculations  were  performed  in  2D  planar  symmetry  on  a  lOOx  100  grid  of  square  15  m  DSMC  cells.  The  DSMC 
time  step  is  4x  1 0  4  s,  the  ratio  of  real  to  simulated  particles  is  1017,  and  the  total  number  of  simulated  particles  was 
approximately  1.3  x  106.  Both  O  and  CO  flow  into  the  cell  at  5656  m/s  from  the  left  and  bottom,  respectively,  at  a 
translational  temperature  of  4  K  and  a  density  of  2x  1016  molecules  per  m3.  The  calculations  were  allowed  4000 
DSMC  steps  to  come  to  steady  state,  after  which  another  4000  steps  were  used  for  sampling.  Each  case  used  8  AMD 
Athlon  CPUs  and  ran  for  roughly  60  minutes.  In  most  cases  (N,  and  A-E)  there  is  no  chemistry.  In  the  remainder  of 
the  cases,  vibrational  transitions  are  treated  as  chemical  reactions  using  the  chemical  mechanism  described  previously. 
True  chemical-change  reactions,  such  as  of  O+O  to  yield  O2  and  O+CO  to  yield  CO2,  were  neglected. 


Influence  of  the  Inelastic  Fraction 


The  first  parameter  examined  is  f ,  the  fraction  of  collisions  that  are  inelastic.  Figures  2  and  3  show  the  translational 
temperature  and  the  number  density,  respectively,  of  the  flow  fields  for  inelastic  fractions  f  of  0.0,  0.1,  and  0.5.  The 
higher  values  of  /,-,  0.5  and  1.0,  have  virtually  identical  results  in  both  temperature  and  density,  suggesting  that  at 
fi= 0.5  the  internal  modes  of  the  molecules  have  already  equilibrated  with  their  local  translational  environment.  At 
fi  of  0,  a  condition  equivalent  to  CO  behaving  as  a  monatomic  gas,  the  flowfield  is  different  than  for  f  >0.5.  The 
most  pronounced  effect  is  seen  in  the  overall  translational  temperature:  the  CO  that  cannot  transfer  energy  into  its 
internal  modes  deposits  it  into  translation,  resulting  in  a  higher-temperature  (14000  K)  contour  along  the  CO  side 
of  the  interaction  zone  than  in  the  other  cases  (10000-12000  K).  The  shape  of  the  density  peak  is  also  affected.  At 
fi  =  0.1  the  results  appear  intermediate  to  those  at  /,=  0.0  and  0.5.  The  composition  of  the  flow  in  all  runs  is  relatively 
constant,  with  a  typical  pattern  shown  in  Figure  1 . 

This  result  is  consistent  with  our  general  expectation  that  the  rotational  temperature  is  often  in  equilibrium  with  its 
local  environment.  It  also  shows  that,  in  this  regime,  the  results  are  not  sensitive  to  the  value  of  f.  On  the  other  hand, 
when  vibrational  degrees  of  freedom  are  lumped  in  with  the  rotational  degrees  of  freedom,  an  inconsistency  could 
develop. 


Influence  of  the  Viscosity  Exponent 

The  second  parameter  examined  is  the  viscosity  exponent  co,  which  was  increased  from  0.75  to  0.90.  Figure  4,  shows 
the  temperature,  number  density,  and  CO  rotational  temperature,  respectively,  in  the  flowfield  at  co  values  of  0.75  and 
0.90.  The  largest  effect  is  in  the  density,  where  the  peak  density  region  becomes  smaller  and  more  symmetrical  for 
the  co=0.90  case.  The  differences  in  terms  of  translational  and  rotational  temperature  are  small.  The  co=0.90  case  has  a 
lower  overall  rotational  temperature,  but  appears  to  have  more  evenly-distributed  translational  temperature. 

What  is  interesting  here  is  that  there  is  a  similarity  in  the  morphology  of  the  interaction  region  for  co=0.90  and  that  at 
co=0.75  with  f  of  0.1  and  0.5.  The  data  suggest  that  an  intermediate  value  of  f  may  have  the  same  effect  on  interaction 
layer  shape  as  a  change  in  co.  Although  variation  in  co  has  been  previously  shown  to  change  the  plume  shape  [12],  the 
effect  of  fj  has  not  been  recently  investigated  in  a  radiation  context. 


Influence  of  the  Internal  DOF 

The  next  examination  compares  CO  frozen  in  the  ground  state  (^=2.0)  with  vibrationally-excitable  CO  with  several 
kinds  of  models:  a  £=4.0  CO  molecule,  a  £=2.0  CO  molecule  with  a  2-state  vibrational  chemistry  model,  and  a 
£=2.0  CO  molecule  with  a  multi-state  vibrational  chemistry  model.  Figures  5  and  6  show  the  number  density  and 
CO  rotational  temperature,  respectively,  for  these  four  cases.  Adding  the  vibrational  DOF  produces  a  denser  CO-O 
interaction  layer,  although  the  flow  for  the  two-state  model  is  not  as  dense  as  that  of  the  multistate  model  and  the  £=4.0 
model.  The  multistate  model  and  the  £=4.0  model  show  densities  that  are  10%  higher  than  for  the  two-state  model, 
and  20%  higher  than  the  vibrationally  frozen  case.  As  one  might  expect  from  the  increased  density,  the  translational 
temperatures  (not  shown)  are  accordingly  lower.  On  the  other  hand,  the  rotational  temperature  of  CO  is  cooler  than 
nominal  in  the  £=4.0  case,  but  much  warmer  in  the  two-state  and  multistate  cases,  with  the  two-state  case  showing  the 
highest  rotational  temperatures. 


Influence  of  Radiation 

To  include  radiation,  the  first-order  decay  reaction  of  CO  at  v  =  1  through  v  =  10  were  included  in  the  calculation. 
We  used  the  Einstein  A  values  for  J  «90  from  Figure  4  of  Braunstein  and  Duff  [26]  for  the  multistate  model  and 
arbitrarily  assigned  a  rough  mean  A  of  60  s~ 1  for  the  2-state  model. 

In  the  previous  section  we  saw  the  different  effects  of  a  2-state  and  a  multistate  model  for  excitation.  The  addition 
of  radiation  does  not  make  a  large  difference  in  this  scenario.  The  general  trend  is  to  reduce  the  flow  field  temperature 
(Figure  7);  this  was  more  prominent  in  the  two-state  model  than  in  the  multistate  model.  There  was  little  observable 
change  in  the  density.  There  was  also  little  observable  change  in  the  CO  rotational  temperature. 


FIGURE  1. 


o 


Simulation 

Domain 


t  f  t  f  f 

CO 


(left)  Schematic  of  simplified  cross-flow  calculation,  (right)  CO  mole  fraction  field  at  /,-  values  of  0.5  (typical). 


FIGURE  2.  Overall  translational  temperature  (K)  fields  at  J)  values  of  0.0,  0.1,  and  0.5. 
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FIGURE  3.  Total  number  density  (m  3)  fields  at  /,-  values  of  0.0,  0.1,  and  0.5. 
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FIGURE  4.  Front  left  to  right:  overall  translational  temperature  (K),  total  number  density  (m  3),  and  CO  rotational  temperature 
(K)  fields  at  ft)  values  of  0.75  and  0.90.  The  top  row  corresponds  to  ft)  of  0.75;  the  bottom  is  ft)  of  0.90. 


FIGURE  5.  Total  number  density  (m  3)  fields.  The  cases  are  for  £=2.0,  £=4.0,  £=2.0  with  a  2-state  vibrational  model,  and  £=2.0 
with  a  multistate  vibrational  model. 


FIGURE  6.  CO  rotational  temperature  (K)  fields.  The  cases  are  for  £=2.0,  £=4.0,  £=2.0  with  a  2-state  vibrational  model,  and 
£=2.0  with  a  multistate  vibrational  model. 


SUMMARY  AND  FUTURE  WORK 

In  this  work  the  sensitivity  of  the  spatial  maps  of  density  and  temperature  are  assessed  in  terms  of  the  VHS  to  parameter, 
the  fraction  of  inelastic  collisions  ft,,  and  the  model  for  the  internal  state  of  the  molecule  through  changes  in  £  or 
through  explicit  chemical  reactions.  Comparing  the  nominal  case  with  the  cases  with  explicit  vibrational  chemistry, 
we  find  that  turning  off  a  reaction  set  (such  as  might  be  done  when  exploring  kinetic  mechanism)  can  significantly 
influence  the  peak  translational  temperature  and  the  rotational  temperature  distribution  in  the  flow  field.  Additionally, 
there  is  evidence  to  suggest  a  covariance  between  the  to  parameter  and  the  in  influencing  the  morphology  of  the 
interaction  layer.  Absent  was  the  a  parameter  in  the  variable  soft  spheres  [27]  potential,  which  may  also  affect  plume 
shape,  and  the  survey  of  different  collision  number  techniques  that  can  effectively  assign  ft  as  a  function  of  local 
conditions,  which  is  relegated  to  future  work. 


FIGURE  7.  Overall  translational  temperature  (K)  fields.  Top  row  shows  nonradiative  results,  bottom  row  shows  radiative-model 
results.  Left  column  is  for  the  2-state  system,  right  column  for  the  multistate  system.  Nominal  Case  N  (no  vibrations)  is  shown  at 
the  far  left  for  reference. 
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